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Abstract. We prove that Dirichlet stereohedra for non-cubic crystallographic 
groups in dimension 3 cannot have more than 80 facets. The bound depends 
on the particular crystallographic group considered and is above 50 only on 9 
of the 97 affinc conjugacy classes of them. 

We also construct Dirichlet stereohedra with 32 and 29 facets for a hexag- 
onal and a tetragonal group, respectively. 



Introduction 

This is the second in a series of three papers (see |2|, devoted to bound the 
number of facets that a Dirichlet stereohedron in 3-space can have. 

A stereohedron in dimension d is any convex polytope which tiles the Euclidean 
space M. d face-to-face under the action of a crystallographic group. A Dirichlet 
stereohedron for a crystallographic group G is the Voronoi region Votgp(P) of a 
point P in the Voronoi diagram of the orbit GP. 

Finding a good upper bound for the number of facets of d-dimcnsional stereo- 
hedra is mentioned as an important problem in |t]] and JlC[ ] , and it is related to 
Hilbert's 18 th problem The most relevant previous results on this problem arc: 

• Delone [Q proved that a stereohedron of dimension d for a certain crystal- 
lographic group G cannot have more than 2 (a + 1) — 2 facets, where a is 
the number of aspects of G, i.e., the order of the quotient group of G by its 
translational subgroup. In dimension 3, crystallographic groups can have up 
to 48 aspects, which gives an upper bound of 390 for the number of facets of 
3-dimensional stereohedra. Non-cubic groups without reflections have up to 
16 aspects (see Table ||), giving a bound of 134. 

• The 3-dimensional stereohedron with the maximum number of facets known 
so far has 38 facets, and it was found by Engel (see (f| and [Q, p. 964]). It is 
a Dirichlet stereohedron for a cubic group with 24 aspects. 

Delone's result contrasts with the fact that arbitrarily large families of non- 
overlapping congruent convex 3-polytopes exist in which the intersection of every 
two polytopes is a common facet ||. 

The main results of this series of papers are summarized in Table |l|. We divide 
the 219 affine conjugacy classes of 3-dimensional crystallographic groups in three 
blocks: Those which contain reflection planes (100 conjugacy classes) were studied 
in d), non-cubic groups without reflections (97 classes) are the object of this paper 
and cubic groups without reflections (22 classes) will be studied in || ■ Our methods 
provide different upper bounds on the number of facets for each group. The two 
columns in Table |l| show respectively the global upper bound obtained for each 
block of groups and the number of affine conjugacy classes of groups where our 
bound is greater than the number of facets of Engel's stereohedron. 

The results for groups with reflections are specially good, since the upper bound 
of 18 facets was proved to be tight in B. For cubic groups we have indicated the 
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global bound 


nbr. of "bad" groups 


Groups with reflections 
1 


18 





Non-cubic groups 
[This paper] 


84 


21 


Cubic groups 


162 


18 


Total 


162 


39 



Table 1. Summary of our results 



results contained in |lj] , which might be slightly improved in the final writing of |3j . 

Table | gives a more detailed description of the results in this paper. It lists the 
21 non-cubic groups where our bound is greater than 38, together with the specific 
bound for each group. The last column in the table indicates where in this paper 
the bound is proved. For some of the groups we have taken Delone's upper bound 
of 8a + 6. 



Group 


Aspects 


Planar group 


Our upper bound 


. . . proved by 


74c2 


8 


P99 


40 


Cor. |3_5| 


ph 1 1 

n g c 


16 


P99 


40 


Cor. pi 


R3 


6 


P 3 


42 


Prop. 


2.7 




P32 


6 


P 3 


42 


Prop. 


2.7 




R3c 


6 


P 3 


42 


Prop. 


2.7 






8 


pgg 


44 


Prop. 


3.4 




P3i2 


6 


P i 


48 


Cor. |L6| 


P3il2 


6 


P i 


48 


Cor. |T| 


P6i 


6 


P i 


48 


Cor. [T| 


P4i22 


8 


P i 


50 


Cor. U 


(J- - - 

a c c 


8 


P 2 


50 


Cor. [T| 


j2 2 2 

a c c 


8 


pgg 


50 


Cor. [T| 


P4i2i2 


8 


P i 


64 


Cor. [T| 


jh 

9 


8 


P 2 


70 


Delone 


P4i22 


8 


P 2 


70 


Delone 


P42d 


8 


P 2 


70 


Delone 


pill 

d d d 


8 


P 2 


70 


Delone 


P6 2 22 


12 


P 2 


78 


Cor. ' 


Ll 


P6i22 


12 


pi 


78 


Prop. 


4.1 




J?3f 


12 




79 


Cor. |2j 


jhl 1 

q c d 


16 


pgg 


80 


Cor. m 



Table 2. Non-cubic groups where our upper bound is more than 38 



Theorem 1 Dirichlet stereohedra for non-cubic crystallographic groups in dimen- 
sion 3 cannot have more than 80 faces. 
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They can not have more than 70 except perhaps in the groups P3§ ; P6i22, P6 2 22 
and Ig-^- They can not have more than 50 except perhaps in these four and five 
other groups. 

The upper bounds are proved in the following Sections, starting from more gen- 
eral methods to more specific ones. For example, Corollary 1.6 gives an upper 
bound derived from simple parameters of each group. This bound is globally 106, 
and it is higher than 38 in only 34 groups. The rest of the paper concentrates in the 
top and bottom parts of Table || lowering the upper bound for the bad groups in 
the bottom of the table and removing groups from the top of the table by showing 
that they cannot produce more than 38 facets. 

Section 4 is essentially devoted to only one group, the group P6i22. We first 
lower the upper bound of 96 from Corollary Ftlj] to a 78, and then show the existence 



of a stereohedron with 32 facets for this group (Example L4). This is the highest 
number of facets of a stereohedron for a non-cubic group obtained to date. 

This example shows, on the one hand, that the general upper bound of 80 is 
at most 2.5 times worse than the actual maximum. And on the other hand, that 
non-cubic groups can produce stereohedra almost as complicated as cubic groups. 
It is worth noting that our construction method shows why it is natural to expect 



stereohedra with "many" facets for the group P6i22. In fact, Lemma 4.2 implies 
that any group whose name (in the International Crystallographic Notation, see ||]) 
contains the string ii2, produces Dirichlet stereohedra with at least 4i + 1 facets. 



We show a stereohedron with 29 facets for the group P4i22 (Example 4.5 ) 



1. A FIRST UPPER BOUND FOR EACH CRYSTALLOGRAPHIC GROUP 

All throughout the paper G will denote a 3-dimensional crystallographic group 
without reflections and P € M 3 will be a base point for an orbit GP, so that the 
Dirichlet stereohedron we want to study is the (closed) Voronoi region Vor<3p(P). 
We assume that P has trivial stabilizer under the action of G, since otherwise GP 
is also an orbit of a proper subgroup of G. Further, there is no loss of generality in 
assuming both P and the metric parameters of G (the parameters which identify 
G in its afhne conjugacy class of groups) to be sufficiently generic since a small 
perturbation to either P or G can only increase the number of facets in the Dirichlet 
stereohedron induced. We say that a point Q S GP is a neighbor of P if Voicp(P) 
and Vor<3p(Q) share a facet. Hence, the number of facets of Vor<3p(P) equals the 
number of neighbors of P. 

1.1. Outline of our method. Let v G M 3 be a vector such that the correspond- 
ing translation t v is in G. We will call vertical the direction of v and horizontal 
the planes orthogonal to it. Let Zp denote the closed infinite band of width 2\v\ 
bounded by the horizontal planes passing through t v P and t_ v P. 

Lemma 1.1. All the neighbors of P lie in Zp. Moreover, the only points in the 
boundary of Zp which can be neighbors of P are t v P and T- V P. 

Proof: Let Q £ GP. Suppose that Q is outside Zp or that it is on its boundary but 
it equals neither t v P nor t^ v P. Without loss of generality assume that PQ ■ v is 
positive (otherwise change v to — v). Consider the points P' = t v P and Q' = t^ v Q. 
The four points PP'QQ', in this order, are the vertices of a planar quadrilateral 
whose angles at P' and Q' are greater or equal than 90 degrees. Since the four 
points are in GP, P and Q cannot be neighbors. □ 
In Zp there are only a finite family ai, . . . ,at of horizontal planes containing 
points of GP. Our method consists on bounding the number of neighbors of P in 
each <Xi separately and using the sum of the numbers obtained as a bound for the 
total number of neighbors. In all cases of interest to us the direction will be chosen 
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so that the subgroup Go of G consisting of horizontal motions (motions which send 
every horizontal plane to itself) contains two independent translations. That is to 
say, it is a 2-dimensional crystallographic group acting on each horizontal plane. 

Since we are assuming P to be sufficiently generic, all the intersections GPHai, 
i € {1, . . . , k} are orbits of Go- The following statement is Lemma 1.3 in 0]. 

Lemma 1.2. Under these assumptions, given two different horizontal planes on 
and aj, a necessary condition for Q € on PI GP = GqQ to be a neighbor of P € 
CKjHGP — GqP is that the two 2-dimensional Dirichlet stereohedra VorQ Q p(P) and 
Vorc Q(Q) overlap (we are implicitly projecting GqP and GqQ along the direction 
of v to a common horizontal plane). □ 

If, moreover, every element of G sends horizontal planes to horizontal planes, 
then we have the following result. We recall that the normalizer of a group Go 
in a bigger group H is the subgroup of H consisting of elements g € H such that 
g~ 1 G g = G . 

Lemma 1.3. Let G be a 3-dimensional crystallographic group. Let a vertical di- 
rection be chosen so that every element of G sends horizontal planes to horizontal 
planes. Let Go be the subgroup of horizontal motions of G. 

Then the orthogonal projection of any orbit ofG to a horizontal plane is contained 
in an orbit of the normalizer of Go in the group of affine isometries of that plane. 

Proof: Let g € G and let h G Go- Since g sends horizontal planes to horizontal 
planes and h sends horizontal planes to themselves, g~ Y hg sends horizontal planes 
to themselves. Hence g~ 1 hg G Go (i.e. Go is normal in G). This, together with 
the fact that g composed with the projection to any horizontal plane is an affine 
isometry on that plane, implies the statement. □ 
The following result, most of which comes from j^], gives a first upper bound on 
the number of neighbors of P in each on, depending on the type of the group Go: 

Theorem 1.4 (|2|, Theorem 3.1]). Let Go be a planar crystallographic group and 
let GqP and GqQ be two orbits of points with trivial stabilizer. Then, the number 
of Dirichlet regions of one of the orbits overlapped by each Dirichlet region of the 
other orbit is bounded above by: 

(i) one, for groups of types pmm, p3ml, pAm and p6m. 

(ii) two, for groups of types cm, cmm, p31m and pAg. 

(iii) four, for groups of types pi, p3, p4, p6, pmg and pm. 

(iv) seven, for the groups p2 and pg. 

(v) eleven, for the group pgg. 

(vi) seven, for the group pgg, if GqP U GqQ is contained in an orbit of the nor- 



malizer of Gq in the affine isometries ofM 2 (see Proposition 3.i ). 
The bounds are all tight except perhaps those for pg and pgg. 

Proof: Parts (i) to (iv) appear in g as Theorem 3.1. Part (v) was also proved in 
0, although not explicitly stated as a result (see Remark 3.2 and the paragraph 



before Lemma 3.5 in that paper). Part (vi) is proved as Proposition 3.3 in Section 
3 of this paper. □ 
We conjecture that the upper bound for pgg is seven even without the normalizer 
condition. 

1.2. The classification of 3-dimensional crystallographic groups. The first 
two invariants to affinely classify a crystallographic group G are its translational 
subgroup and its point group. The translational subgroup T is the subgroup con- 
sisting of translations, and belongs to one of the 14 Bravais lattice types. The 
point group of G is the quotient group [G : T}. The elements of [G : T] are the 



NUMBER OF FACETS OF THREE-DIMENSIONAL DIRICHLET STEREOHEDRA II 5 



aspects of G. The point group can be understood as a discrete group of motions 
in the 2-dimensional unit sphere satisfying the "crystallographic restriction" : there 
are no rotations of order 5 or greater than 6. There are 32 groups satisfying this, 
modulo isometries of the sphere, and all of them appear as point groups of some 
crystallographic group. 

In turn, the lattices are classified into seven crystallographic systems: monoclinic, 
triclinic, orthorhombic, hexagonal, trigonal, tetragonal and cubic. This classification 
is according to the point groups of their normalizers in the group of affine isometries 
of K 3 , except that most authors include in the trigonal system some groups whose 
lattice is hexagonal, the distinction between the hexagonal and trigonal systems 
relying on properties of the point group. 

Groups which contain reflection planes were studied in Q and cubic groups 
will be dealt with in jj). Table || lists the 97 affine conjugacy classes of non- 
cubic 3-dimensional crystallographic groups without reflections, divided according 
to crystallographic systems and point groups. We use the International Crystallo- 
graphic Notation to name 2-dimensional and 3-dimensional crystallographic groups 
and have taken the monograph || as a source for their classification. The capital 
letter at the beginning of a name indicates the Bravais type of the lattice within 
its crystallographic system. The rest of the name encodes generators for the group 
in a certain way. 

Only 58 of the 97 groups in Table || have more than 4 aspects. We will deal only 
with these 58 groups since the Delone bound for groups with four or less aspects is 
already 38, matched by Engel's example. Moreover, experimental evidence indicates 
that in groups with few aspects Delone's bound is a good approximation to what 
can happen in the worst case: There are Dirichlet stereohedra for groups with one 
and two aspects having 14 and 20 facets respectively (?], p. 963]. 

The 58 groups belong to the following crystallographic systems: orthorhombic 
(10 groups), trigonal (13 groups), hexagonal (9 groups) and tetragonal (26 groups). 
They are listed in Tables 

There are the following 8 possible Bravais types of lattices in the 58 groups of 
Tables g-[jj 

• The three primitive lattices (noted P in Table ||), generated by a vertical 
translation and the following planar lattices in horizontal planes: a rectangu- 
lar lattice in the orthorhombic system, a triangular lattice in the hexagonal 
and trigonal systems, and a square lattice in the tetragonal system (see Fig- 
ure |l|). Each primitive lattice has a primitive cell associated to it, which is 
an orthogonal prism over a rectangle, triangle or square respectively. The 
primitive cell tiles the space face-to-face producing as vertices of the tiling an 
orbit of the lattice. 



• • • 

• • • 



• ••• 

rectangular square triangular 

Figure 1. The three planar lattices 

The two body-centered lattices (noted /) in the orthorhombic and tetragonal 
systems. They are generated by the primitive lattice and any translation 
sending a vertex of the primitive cell to the centroid of the cell. Hence, they 
contain the primitive lattice as a sublattice of index 2. 
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System 


Point (aspects) 
group 


Crystallographic groups 
without reflections 


Triclinic 


1 (1) 


PI 


1 (2) 


PT 


Monoclinic 


2 (2) 


P2 P2i B2 


m (2) 


Pb Bb 


- ( 4 ) 


P2 p2i S| 

f» f» b 


Orthorhombic 


222 (4) 


P22 P2i22 P22 x 2i P2i2i2i 
C222 C2i22 F222 J222 7222' 


2mm (4) 


P2cc P2ica P2cn P26a P2ifen 
P2nn C2cc /2cc F2dd A2ba 


2.2.2. (g) 


pAAA pAA±x P i zi ±1 

n n n a n n n c c 

P A AL i. P AL AL AL P 1 _ £l 

acc a & c nca 

ioo , O 1 o _ooo ooo 

r> Z Z z AAA T A A A rp A A A 

n b a acc acc d d d 


Trigonal 


3 (3) 


P3 P3i P3 


3 (6) 


P3 P3 


32 (6) 


P32 P3i2 P312 P3il2 P32 


3m (6) 


P3c P3ic P3c 


3^ (12) 


P32 P3l2 P32 

C C C 


Hexagonal 


6 (6) 


P6 P6i P6 2 P63 


622 (12) 


P622 P6i22 P6 2 22 P6 3 22 


6mm (12) 


P6cc 


Tetragonal 


4 (4) 


P4 P4i P4 2 14 P4i 


4 (4) 


P4 J4 


£ (8) 


p4 pAz. J a\i 

n n q 


422 (8) 


P422 P4i22 P4 2 22 P42j2 
P4i2i2 P4 2 2i2 J422 74i22 


4mm (8) 


P4 2 gc P4nc P4cc 74icd 


42m (8) 


P42c P42ic P4g2 
P4c2 P4n2 P4c2 J42d 


(16) 


p£2i2 pla.22 p£21 j4i 2 2 

n c c n q c n n c q c d 



Table 3. The 97 non-cubic 3-dimensional crystallographic groups 
without reflections 



The face-centered lattice in the orthorhombic system (noted F), generated 
by the primitive lattice and translations from a vertex of the primitive cell to 
the centroids of the three facets of the cell incident to that vertex. It contains 
the primitive lattice as a sublattice of index 4. 

The base-centered lattice in the orthorhombic system (noted C), generated 
by the primitive lattice and the translation from a vertex of the primitive cell 
to the centroid of one of the three facets of the cell incident to that vertex. It 
contains the primitive lattice as a sublattice of index 2. 

The rhombohedral lattice in the trigonal system (noted R), generated by the 
primitive lattice and a translation from one of the vertices of the primitive 
cell (a triangular prism) to a point in the axis of the primitive cell and whose 
distances to the two triangular bases of the prism are in the ratio 1:2. It 
contains the primitive lattice as a sublattice of index 4. 
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Group 


a 


Go 


a 


i 


Z 


Cor. L6| 


Final bound 


Remarks 




pill 

n n n 


8 


p2 


2 


7 


1 


50 


38 


Prop. 


2.5 




plllx 

a n n 


8 


P 2 


2 


7 


1 


50 


38 


Prop. 


2.5 


p2 2\ 2i 

n c c 


8 


P 2 


2 


7 


1 


50 


38 


Prop. 


2.5 


p2 2i 2 

a c c 


8 


P 2 


2 


7 


1 


50 


38 


Prop. 


2.5 


p2±2±2± 

a b c 


8 


pg 


2 


7 


1 


50 


38 


Prop. 


3.7 




p2j_ 2 2i 

n c a 


8 


pg 


2 


7 


1 


50 


38 


Prop. 


3.7 




pill 

n b a 


8 


pgg 


4 


7 


1 


22 










(Jill 

a c c 


8 


P 2 


2 


7 


1 


50 










Jill 

a c c 


8 


pgg 


4 


7 


2 


50 










pill 

d d d 


8 


P 2 


2 


7 


2 


106 


70 


Delone 





Table 4. Orthorhombic groups without reflexions with more than 
4 aspects 



Group 


a 


Go 


a 


i 


2 


Cor. f| 


Final bound 


Remarks 


P3 


6 


P 3 


3 


4 


1 


16 








RS 


6 


P 3 


3 


4 


1 


48 


42 


Prop. 


2.7 


P32 


6 


pi 


3 


4 


1 


16 








P3i2 


6 


pi 


1 


4 


1 


48 








P312 


6 


P 3 


3 


4 


1 


16 








P3il2 


6 


pi 


1 


4 


1 


48 








P32 


6 


P 3 


3 


4 


3 


48 


42 


Prop. 


2.7 


P3c 


6 


P 3 


3 


4 


1 


16 








P31c 


6 


P 3 


3 


4 


1 


16 








i?3c 


6 


P 3 


3 


4 


3 


48 


42 


Prop. 


M 


P3^ 

C 


12 


P 3 


3 


4 


1 


32 








P31§ 


12 


P 3 


3 


4 


1 


32 








m( 


12 


p3 


3 


4 


3 


96 


79 


Cor. 


2.8 



Table 5. Trigonal groups without reflexions with more than 4 aspects 



Group 


a 


Go 


a 


i 


2 


Cor. i.e 


Final bound 


Remarks 


P6 


6 


p6 


6 


4 


1 


8 








P6i 


6 


pi 


1 


4 


1 


48 








P6 2 


6 


p2 


2 


7 


1 


36 








P6 3 


6 


P 3 


3 


4 


1 


16 








P6cc 


12 


p6 


6 


4 


1 


16 




P622 


12 


p6 


6 


4 


1 


16 








P6 3 22 


12 


P 3 


3 


4 


1 


32 








P6i22 


12 


pi 


1 


4 


1 


96 


78 


Prop. 


4.1 


P6 2 22 


12 


p2 


2 


7 


1 


78 









Table 6. Hexagonal groups without reflexions with more than 4 aspects 
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Group 


a 


Go 


a 


i 


Z 


Cor. 1.6 


Final bound 


Remarks 




n 


8 


pA 


A 


A 


1 


16 


— 








ph. 

n 


8 


P 2 


2 


7 


1 


50 


38 


Prop. 


2.5 


a 


8 


P 2 


2 


7 


1 


106 


70 


Delone 




P422 


8 


pA 


A 


4 


1 


16 


— 








P4 X 22 


8 


pl 


1 


4 


1 


64 


50 


Cor. f(\ 




P4 2 22 


8 


p2 


2 


7 


1 


50 


38 


Prop. 


2.5 


P42i2 


8 


pA 


4 


4 


1 


16 


— 








P4i2i2 


8 


pl 


1 


4 


1 


64 


— 








P4 2 2i2 


8 


P 2 


2 


7 


1 


50 


38 


Prop. 


2.5 


7422 


8 


pA 


4 


4 


2 


32 


— 








P4i22 


8 


p2 


2 


7 


2 


106 


70 


Delone 




P4 25 c 


8 


P99 


4 


7 


1 


22 


— 








P4tt,c 


8 


pA 


4 


4 


1 


16 


— 








P4cc 


8 


pA 


4 


4 


1 


16 


— 








M-Lcd 


8 


P99 


4 


7 


2 


50 


44 


Prop. 


3.4 




PA2c 


8 


P 2 


2 


7 


1 


50 


38 


Prop. 


2.5 




P42ic 


8 


P 2 


2 


7 


1 


50 


38 


Prop. 


2.5 


P4. 9 2 


8 


P99 


4 


7 


1 


22 


— 








P4c2 


8 


P 2 


2 


7 


1 


50 


38 


Prop. 


2.5 


PAn2 


8 


P 2 


2 


7 


1 


50 


38 


Prop. 


2.5 


IAc2 


8 


P99 


4 


7 


2 


50 


40 


Cor. U 




IA2d 


8 


P 2 


2 


7 


2 


106 


70 


Delone 




p± 2i 2 

nee 


16 


pA 


4 


4 


1 


32 










p422 

n n c 


16 


pA 


4 


4 


1 


32 










pig. 2 2 

n g c 


16 


P99 


4 


7 


1 


50 


40 


Cor. [ 






q c d 


16 


P99 


4 


7 


2 


106 


80 


Cor. [ 







Table 7. Tetragonal groups without reflexions with more than 4 aspects 



1.3. A first upper bound. Suppose that G is one of the 58 groups of Tables | @. 
The lattice of G contains as a sublattice one of the three primitive lattices, and we 
choose as vertical the direction of the principal axis in the corresponding primitive 
cell. The principal axis is well defined in the hexagonal, trigonal and tetragonal 
systems. In the orthorhombic system there are three perpendicular and equivalent 
axes. We choose one arbitrarily unless the lattice is base-centered (noted C) in 
which case we choose as vertical the direction perpendicular to the centered faces. 

Go is crystal- 
We 
The 



This choice fulfills the conditions required in Lemmas 1.2 and 1.3 



lographic and every element of G sends horizontal planes to horizontal planes, 
take as vector v to apply Lemma 1.1 the shortest vertical translation in G. 
following Lemma shows how to compute the number of horizontal planes in the 
band Zp of Lemma 1.1 which contain points of GP. This number depends only on 



the type of the lattice and the number of aspects of G and Go- 

Lemma 1.5. Apart of the two boundary planes and the plane containing P , there 
are exactly 2al/a^ — 2 horizontal planes in Zp containing points of the orbit GP, 
where: 
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• a is the number of aspects of G . 

• ao is the number of aspects of the horizontal subgroup Gq . 

• I is a number depending only on the lattice of G, and takes the value 1 in prim- 
itive and base-centered lattices, 2 in body- centered and face- centered lattices, 
and 3 in rhombohedral lattices. 

Proof: Let T and T denote the translational subgroups of G and Go, respectively. 
The number I is the ratio of horizontal planes containing points of T to planes 
containing points of the primitive lattice, so that each orbit of T projects to I orbits 
ofT . 

The closed band Zp contains 2al plus ao orbits of To- Each horizontal plane 
containing points of GP contains ao orbits of To. Hence, there are 2al/ao + 1 such 
planes, including the two boundary planes and the one containing P. Subtracting 
3 gives the statement. □ 

Corollary 1.6. Let i be the maximum number of intersections between two differ- 



ent orbits of the planar crystallographic group Go, as stated in Theorem 1.4- Then, 
the number of neighbors of the point P in an orbit GP is bounded above by 

2i{al/a - 1) + 8, 

where a, ao and I are as in Lemma 1 



Proof: By Theorem 1.4, P has at most i neighbors in each horizontal plane of Zp 
other than the two boundary ones and the one containing P. The only possible 
neighbors in the boundary are t v (P) and T— V (P), and in the plane containing P 
there are at most six neighbors, since a planar Dirichlet stereohedron has at most 
six facets. □ 



In Tables [|-[^, the column "Cor. L6" indicates the upper bound given by this 



statement in each of the 58 groups of interest to us. 

2. Influence region. Groups with a horizontal p2 or p3 



Corollary |l.q is based on using global bounds for the number of regions of Vorc Q 
which can intersect the region Vor<3 p(P). But of course the number of regions in- 
tersected depends on Q, and we are interested only in a few concrete possibilities for 
the point Q in each 3-dimensional crystallographic group. To study the dependence 
of the number of intersections with Q we recall the formalism of extended Voronoi 
regions and influence regions, already used in Q . 

2.1. Influence regions in planar crystallographic groups. Let Go be a crys- 
tallographic group in the plane. Let A<o be a crystallographic group containing Go- 
We will call fundamental subdomains of Go (w.r.t. No) the fundamental domains of 
./Vo- For a fixed fundamental subdomain D of Go, we call extended Dirichlet region 
of D under Go (and denote it Extc (-D)) any region containing UQ e _DVorG Q(Q). 

Observe that, for any given D, the fundamental subdomains {tD : t £ Nq} tile 
the plane. We call influence region of D under Go (and denote it lnfio (D)) the 
union of all the tiles tD with r G No for which ExtG (rD) and Extc (D) overlap. 

Lemma 2.1. Let P and Q be two points in the plane and let D be a fundamental 
subdomain containing P. A necessary condition for Vorc p{P) and Vorc Q{Q) to 
overlap is that Q lies in Infl G (D). 

Proof: Let t e N be such that Q S tD. By construction, Votg p(P) C ExtG (-D) 
and VorG Q(Q) C ExtG (TD). This proves the statement. □ 
Observe finally that the action of Go divides the fundamental subdomains in a 
finite number of orbits. Then, for each point Q € K 2 , the number of regions of 
Vor<3 Q which intersect Votg p{P) can be bounded above by counting how many 
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fundamental subdomains in the orbit of fundamental subdomains containing Q 
intersect the influence region InfiG (Z)). 

In this section we apply this formalism to 2-dimensional crystallographic groups 
of types p2 and p3: 

Groups of type p3: Let Go be a group of type p3, generated by an order 3 rotation 
and two translations of equal length with angle ir/6 to one another. The order- 
3 rotation centers of the group form a triangular lattice, and any two adjacent 
elementary triangles of the lattice form a fundamental domain. Our choice of Nq is 
generated by Go together with any reflection on one of the sides of an elementary 
triangle. Each elementary triangle is a fundamental subdomain, and there are two 
Go-orbits of fundamental subdomains, which we call B and W for "black" and 
"white" . The three neighbors of a black triangle are white and vice- versa. We 
assume the subdomain containing P to be white. 

The extended Dirichlet region of an elementary triangle consists of this triangle 
and the three adjacent to it. This is shown in the left part of Figure || in the 
following way: a base point P in a fundamental subdomain D and other six points 
of the orbit Go-P are drawn. Let tP be one of the six points. The region of the 
plane consisting of points closer to tP than to P for any choice of P in the central 
fundamental subdomain is clearly disjoint with ExtG Q (Z3). What is left is precisely 
the union of D and its three adjacent subdomains. 

Hence, the influence region consists of the fundamental subdomains which are 
either in ExtG (_D) or adjacent to a fundamental subdomain in Ext<3 (.D). This 
gives the 10 fundamental subdomains in the right part of Figure |[ 

In total, there are seven white and three black triangles in the influence region. 



/ WW w\ 

*\ B /£\b /N 

w\ / w\ / w 



\ B / 
w\/ w 



Figure 2. The influence region for a group of type p3 




Lemma 2.2. Let Gq be a group of type p3. Let P and Q be two points in the plane 
and suppose that P lies in a white triangle, in the sense explained above. Then, the 
number of regions of Vorc Q which overlap VorGoP(P) is at most three if Q is in 
a black triangle and at most four if it is in a white triangle. 

Proof: The bound of three for black triangles is obvious since there are only three 
black triangles in the influence region. For the white triangles we would in principle 



have a bound of seven, but part (iv) of Theorem 1.4 gives four. □ 



Groups of type p2 with rectangular grid: Let Go be a group of type p2, generated by 
two independent translations and a rotation of order two. We assume further that 
the two generating translations are orthogonal to one another, which implies that 
the order-2 rotation centers form a rectangular grid. Any two adjacent elementary 
rectangles in this grid form a fundamental domain for Go- Our choice of Nq is 
generated by Go together with any reflection on one of the sides of an elementary 
rectangle. Then, each elementary rectangle is a fundamental subdomain and there 
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are two Go-orbits of fundamental subdomains which we call B and W for "black" 
and "white" . The four neighbors of a black rectangle are white and vice- versa. We 
assume the rectangle containing P to be white. 

The extended Dirichlet region of an elementary rectangle consists of this rec- 
tangle and the four adjacent to it. The influence region consists of the extended 
Dirichlet region and the eight elementary rectangles adjacent to it. See Figure 
In total, there are nine white and four black rectangles in the influence region. 
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Figure 3. The influence region for a group of type p2 



Lemma 2.3. Let Go be a group of type p2 with a rectangular grid. Let P and Q 
be two points in the plane and suppose that P lies in a white rectangle, with the 
meaning explained above. Then, the number of regions of Vorc Q which overlap 
VorGoP(P) is at most four if Q is in a black rectangle and at most seven if it is in 
a white rectangle. 

Proof: The bound of four for black rectangles comes from the influence region. The 
bound of seven for white rectangles comes from Theorem 1.4. □ 
We are interested in the case where Go is the horizontal subgroup of a 3- 
dimensional crystallographic group G. We will always make sure that N is a 
group of horizontal motions that not only contains Go but is also contained in the 
normalizer of G in the isometry group of R 3 . The following lemma implies that 
if this happens then there is no loss of generality in assuming that our base point 
P lies over any fixed fundamental subdomain D, since if this is not the case then 
there will be an orbit GP' isometric to GP (in particular, with congruent Dirichlet 
stereohedra) and with P' over D. 

Lemma 2.4. Let G be a 3-dimensional crystallographic group. Let Go be its hor- 
izontal subgroup and let Nq be a horizontal group containing Gq and contained in 
the normalizer ofG. Let D be a fundamental subdomain o/Go- Let P be any point 
in R 3 . Then, there is an isometry t sending P to a point P' over R and with 
t(GP)=G(tP). 

Proof: Let r £ Nq be a horizontal motion sending P to a point P' over D. This 
exists since D is a fundamental domain of Nq. Since r is in the normalizer of G, 
tGP = GtP = GP'. This finishes the proof. □ 

2.2. Groups with a horizontal p2 of rectangular type. 

Proposition 2.5. The four orthorhombic and seven tetragonal groups with 8 as- 
pects displayed in Figure Q have a horizontal group Gq of type p2 with rectangular 
grid. The reflections on vertical planes containing edges of elementary rectangles 
of the grid lie in the normalizer of G. 

In each case, the band Zp contains six horizontal planes (not counting the bound- 
ary and middle ones) of a generic orbit GP, four which produce black orbits and 
two which produce white orbits. Hence, their Dirichlet stereohedra cannot have more 
than 38 facets. 
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Figure 4. Orthorhombic and tetragonal groups whose horizontal 
group is a rectangular p2 



Proof: The bound of 38 follows from the rest of the statement by a counting argu- 
ment, using Lemma |2~3| : to the 6 neighbors in the plane of P and 2 in the boundary 
of Zp we have to add four in each of the four "black" planes of orbits and seven in 
the two "white" ones. 

The rest of the statement can be easily checked in the graphical representation 
of each group displayed in Figure ^. Let us first explain this representation. 

In each picture, a horizontal projection of a generic orbit appears. The black 
ovals in the corners of the rectangles/squares represent the order-2 rotation centers 
of the horizontal group p2. Hence, any of the four rectangles is a fundamental 
subdomain, any adjacent two are a fundamental domain of Go, and the four together 
are a fundamental domain of the translational subgroup of Go- In each picture, the 
black circle in the bottom-right quadrant represents the base point P for the 3-d 
orbit GP. The other seven circles represent points in the other seven aspects of G, 
with their heights indicated in the following way: let z be the height of the base 
point and assume that the shortest vertical translation in G has length 1. Then, 
black circles with a number a represent horizontal orbits at heights z + a (and 
hence at z + a + m for any integer m) while white circles with a number a represent 
orbits at height — z + a (and hence at —z + a + jn for any integer m). The absence 
of a number means a = (as happens in the base point itself). 

To make all this clearer, let us describe generators for the group P— — — deduced 
from its graphical representation. First, there is the translational subgroup, gen- 
erated by horizontal translations on the sides of the big rectangle and a vertical 
translation of length 1. Then, we have to describe how to get the other seven points 
in the figure from the base point. The white point in the bottom-right rectangle 
is obtained by a rotation of order two around the vertical axis of the bottom-right 
rectangle, followed by a reflection on the horizontal plane at height zero. The two 
points on the upper-left rectangle are obtained from the two in the bottom-right 
rectangle by a rotation of order two on the vertical line through the middle black 
oval. The other four points are obtained from the first four by a rotation of order 
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two around the line at height 1 /4 over any of the two horizontal axes of the big 
rectangle. 

Let us also see how to check in the picture the conditions of the statement. 
An isometry r is in the normalizer of G if and only if for any generic orbit GP, 
tGP is again an orbit of G (if this happens, let P' = tP and let P" be such 
that tGP = P". Then tGt^P' = GP" which implies P 1 € GP" and, since P is 
generic, tGt -1 — G). In each figure it is easy to check that a reflection in any of 
the displayed lines produces a new orbit of G. 

As for the number of planes in the band Zp, each horizontal orbit of G, except 
for the one containing P, produces exactly two planes in Zp, since Zp has width 
2 and the minimal vertical translation has length one. In the Figures, the left-top 
and right-bottom quadrants are "white" fundamental subdomains and the other 
two are "black". □ 

More or less the same technique can be applied to the group P4i22, depicted 
in Figure |. Part (a) is its standard graphical representation, showing the eight 
points of a generic orbit GP which lie inside the translational primitive cell of 
the tetragonal system, projected onto the XY plane and with their heights (Z- 
coordinates) indicated. Part (b) shows the projection of these same points to the 
YZ coordinate plane. The black ovals indicate order-2 rotations of G with axes in 
lines parallel to the X axis. They can also be read as order-2 rotation centers in 
the planar subgroup of G which preserves planes parallel to YZ. This subgroup, 
hence, is of type p2 with a rectangular grid. 




Figure 5. The group P4i22 



Corollary 2.6. Dirichlet stereohedra for the group PA\22 cannot have more than 
50 facets. 

Proof: Let G be a group of type PA\22. Let Gi be the subgroup of G which 
preserves planes parallel to YZ. Let Zp be the band centered at the YZ-plane 
containing our base point P and of width two times the shortest translation of G 
(and of Gi) in the X-direction. Only seven Gi-orbits intersect the interior of the 
band. Although we cannot say h ere t hat the Gj-o rbits are relat ed b y the normalizer 
of Gi , still we can apply Lemma [O], Lemma |l.2j and Theorem L4 to conclude that 



each Gi orbit produces at most 7 neighbors except G\P which produces at most 6. 
This gives at most 48 neighbors in the interior of Zp and we have to add the two 
translates of G\ in the ^-direction. 

Observe that the choice of P shown in Figure || produces four Gi-orbits in 
"black" fundamental subdomains, for which the number of neighbors is at most 4. 
But there are choices of P in which all the Gi-orbits lie in white subdomains. □ 



2.3. Groups with a horizontal p3. 
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Proposition 2.7. The four trigonal groups displayed in Figure \Q have a horizon- 
tal group of type p3. The first three have ten horizontal planes (not counting the 
boundary and middle ones) in the region Zp, six of which lie over black subdomains 
and four over white subdomains, in the sense of Lemma 2.i. Hence, their Dirichlet 
stereohedra cannot have more than 42 facets. 

The last group has 22 horizontal planes (not counting the boundary and middle 
ones) in the region Zp, 12 which produce black orbits and ten which produce white 
orbits. Hence, their Dirichlet stereohedra cannot have more than 84 facets. 




Figure 6. Trigonal groups whose horizontal group is of type p3 



Proof: As in Proposition |2.5| , the stated bound follows from the rest of the state- 
ment. The statement on horizontal planes follows from the graphical representation 
of each group in Figure 0. For each group, the small black triangles represent cen- 
ters of rotation of order three. The region displayed is a fundamental domain for 
the horizontal group and the six (in the first three groups) or twelve (in the last 
one) black or white points rep rese nt a generic orbit of G, with the same conventions 
as in the proof of Proposition 2J3. □ 
Fo r the trigonal group i?3| the upper bound of 84 neighbors obtained in Propo- 
sition 2.7 can still be lowered a bit if we take into account that not only horizontal 



planes but also vertical planes contain many points of each orbit. Observe that 
the bound of Proposition 2/7 can be stated more precisely as saying that P can 
have at most 48 neighbors which project to white triangles in the influence region 
and at most 36 neighbors which project to black triangles. We will be interested 
in the latter ones. Figure shows the projection of a generic orbit to the union of 
4 fundamental subdomains forming a triangle. If we assume P to project to the 
central triangle, colored white, then the other three triangles are the three black 
triangles contained in the influence region. 

The dashed lines in Figure represent nine vertical planes A, B, C, A' , B' , 
C , A" , B" and C" with the property that any black neighbor of P lies in two of 
them. On the other hand, the subgroup of G fixing each of those planes is a planar 
crystallographic group of type p2 and, by Theorem 1.4, produces at most seven 
neighbors. Hence the number of black neighbors cannot exceed 9 x 7/2 = 31.5. 
This implies the following: 
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Figure 7. Vertical planes in the group i?3| 

Corollary 2.8. A Dirichlet stereohedron for a group of type i?3^ cannot have more 
than 31+48=79 neighbors. 

3. Reduced influence region. Groups with a horizontal pgg or pg 

3.1. How two orbits of pgg related by the normalizer intersect. We recall 
the following result from [|| (Theorem 3.1): 

Lemma 3.1. Let Go be a planar crystallographic group of type pgg and let P and 
Q be two points in the plane. Assuming that GqP U GoQ is an orbit of a certain 
crystallographic group, the number of Dirichlet regions of Vorc a Q which overlap 
VorG p(P) * s a t most seven. □ 

A case study shows that under the hypotheses of the lemma we must have 
Q = tP where r is an element of the normalizer of Go in the group of Euclidean 
isometries of the plane. We intend to extend the lemma to this more general case, 



thus proving part (vi) of Theorem 1.4 



Recall that pgg is generated by two perpendicular translations 
together with a rotation of order two and a glide reflection of vector half of one 
of the translations. There are two possibilities for the normalizer: 

• If the generating translations have different lengths, then the normalizer iVo 
is generated by Go together with reflections on the lines supporting the rect- 
angles of the grid of rotation centers. Go has index four in iVo and we have 
that for any r € -/Vo, t 2 is in Gq. In particular, Go U rGo is a crystallographic 



group for every r G No, and hence Lemma 3d already implies what we want 
to prove. 

• If the generating translations have the same length, then the normalizer iVo is 
generated by Go together with the reflections mentioned above and rotations 
of order four in the rotation centers of Go- This is the case we will be interested 
in. 

Go has index eight in No. More precisely, we can take as fundamental domains 
of Go squares with vertices in rotation centers and as fundamental subdomains 
(i.e. fundamental domains of Nq) the eight triangles in which the symmetries of 
the square divide the fundamental domains. The influence region is computed in 
Figure @. The left part is an extended Dirichlet region of the initial fundamental 
subdomain D, computed as in Figure ^. The right part is the union of all the 
fundamental subdomains whose extended Dirichlet region overlaps the one on the 
left part. The extended Dirichlet region of a fundamental subdomain tD is obtained 
applying r to the extended Dirichlet region ExtG (-D), for each t E Nq. 
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Figure 8. Influence region of a pgg with a square grid of rotation centers 

We have labeled the eight elements of Nq/Gq (and hence the fundamental sub- 
domains) with the letters A, B, C, D, A', B', C and D'. A represents Go itself and 
A, B, C and D would form the normalizer of Go in case that the two generating 
translations had different length. The number of fundamental subdomains of each 
type in the influence region is shown on the right part of Figure |8|. This number 
never exceeds nine. This implies: 

Lemma 3.2. Let Go be a planar crystallographic group of type pgg whose gener- 
ating translations have equal length. Let GqP and GqQ be two orbits with trivial 
stabilizer. Then, the number of Dirichlet regions of one of the orbits overlapped by 
each Dirichlet region of the other orbit is at most nine. 

So far we have not used the fact that we are interested in orbits related by the 
normalizer. How to use this property is exhibited in Figure ^|. In part (a), our base 
orbit of Go is shown (in black) together with the corresponding orbit of type "B" 
(in white). One point of this latter orbit has been crossed out, meaning that its 
corresponding Dirichlet region can never overlap the region of the base point. The 
reason is in the four points joined by a dashed quadrilateral: the movement of the 
normalizer sending the black vertices of this quadrilateral to the white ones is an 
order-2 rotation, hence the quadrilateral is a parallelogram. The bisectors of the 
two black vertices and the two white vertices of the quadrilateral will be parallel 
and will separate the Dirichlet regions of the base point P and the point crossed 
out. 

The same argument applies to the point crossed out on the top end of Figure 
||(b), where the "A" and "C" orbits are shown. But we have crossed out also two 
other points, one on the top-right end and one on the bottom-left end. The proof 
that the Dirichlet regions of these two points cannot overlap the one of our base 
point is in Q (paragraph [A, C] and, in particular, Figure 19). Hence, instead of 7 
and 9 fundamental subdomains of types B and C which potentially could produce 
Dirichlet regions overlapping the base Dirichlet region we have now 6 of each type. 
As a conclusion: 

Proposition 3.3. Let Go be a planar crystallographic group of type pgg and let P 
and Q be any two points in the plane. Assuming that Q = tP for some element r 
in the normalizer o/Go, the number of Dirichlet regions of Vorc Q which overlap 
VorG p(P) is at most seven. □ 

3.2. Groups with a horizontal pgg. The reduced influence region of Figure 



can be used to lower the upper bounds given by Corollary 1.6 for the groups whose 



horizontal subgroup Go is a pgg with a square grid. This follows from the fact 



NUMBER OF FACETS OF THREE-DIMENSIONAL DIRICHLET STEREOHEDRA II 17 




(c) 

Figure 9. Some fundamental subdomains in the influence region 
in Figure || cannot produce overlapping Dirichlet regions, if the 
two orbits are related by the normalizer (parts (a) and (b)). This 
produces a "reduced influence region" (part (c)). 



that only one type of planar orbits of those produced by the normalizer give seven 
neighbors (the one we have labeled D), and the rest only six. 

More precisely, Figure [lO] shows the graphical representation of the four groups 
in question. The square displayed is a fundamental domain of Go, divided into 
eight fundamental subdomains. In the first three groups the "bad" planar orbit, 
labeled with a D in Figure || does not appear. Hence we can take i — 6 in the 



computations of Corollary |l.6| for these groups. In I^cl ^ ne or t>it appears, 



but still we can count six neighbors for 12 of the planes in the band Zp and seven 
for only two of them. This gives: 




lA x cd 




14 c2 




.4_ 2 22 
n gc 




Figure 10. Tetragonal groups whose horizontal group is of type pgg 
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Proposition 3.4. Dirichlet stereohedra for the groups I A\cd, 74c2 and -P^|f can- 
not have more than 44 facets. Dirichlet stereohedra for the group I— 1 1 cannot have 
more than 94 facets. 

But we can use vertical planes to refine this result a bit. Figure |ll| shows the 
projection of a generic orbit GP to the influence region computed above, where 
G is of type 74c2 or P^-\ in P ar t (a) of the Figure and of type in part 



(b). The number 44 or 94 stated in Proposition 3.4 is obtained counting 2 possible 
neighbors over each fundamental subdomain, except in the six grey subdomains 
(those containing points at the same height as P except the one containing P 
itself) where only one neighbor over each fundamental subdomain is possible. The 
dashed line in part (a) of the figure represents a vertical plane containing the base 
point P where we have counted ten possible neighbors while there are at most six 
real neighbors (since in a planar Dirichlet tiling each region is a neighbor of at most 
other six). Hence, instead of 44 we can take 40 as an upper bound for 74c2 or 

p42.2 2 





(a) 14 c 2 and 



n gc 



(b) 



Figure 1 1 . Vertical planes in the groups of type pgg 

Similarly, the five dashed lines in part (b) of Figure |ll| represent five vertical 
planes A, B, C, D and E where we have counted 12, 12, 10, 10 and 12 possible 
neighbors respectively, giving a total of 50 because two of them are common to A, 
B and C, and other two common to B and D. As before, A, B and C can produce 
at most six neighbors each, because they all contain P. 

Let Gi denote the subgroup of G preserving planes parallel to A, in particular 
preserving D and E. G\ is of type pgg (although it does not have its two translations 
of equal length). According to Theorem 1.4 each plane parallel to A can contain at 
most 11 neighbors, and only seven if the Gi-orbits in that plane and in A are related 
by the normalizer. The latter happens in the plane -D; the orbits in the planes A 
and D are related by a glide reflection in the bisecting plane, which projected to 
any of them becomes a translation lying in the normalizer of the corresponding pgg. 

Hence, the planes A, B, C, D and E can contain in total at most 6+6+6+11+7 = 
36 neighbors, instead of the 50 which we have counted. 

Corollary 3.5. Dirichlet stereohedra for the groups IAc2 and P^ ~ cannot have 
more than 40 facets. Those for -f-^r§§ cannot have more than 80 facets. 
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Remark 3.6. Using the methods of H, Section 3.3] i n a m ore sophisticated way, 



the number appearing for pgg in part (v) of Theorem 1.4 can be lowered to nine 
instead of eleven (see B Remark 3.2]). In the preceeding argument, this would 
lower the number of possible neighbors in the plane E by two as well, giving a 
bound of 78 instead of 80. But this small improvement does not seem to be worth 
the effort of giving a proof here. 

3.3. Groups with a horizontal pg. We will now compute the influence region 
of a planar group of type pg, in order to lower the bound for the groups P^ 1 !^ 1 
and -P^^^r- Remember that pg is generated by two perpendicular translations 
and a glide reflection on a line parallel to one of them. Let a and b be the vectors 
of the two translations, a being parallel to the glide-reflection line. Observe that G 
consists of translations and glide-reflections on a family of lines parallel to a and 
distant V to one another. Any rectangle of sides a and ^ placed between two 
consecutive such lines is a fundamental domain for Go . See Figure ^(a) . 

We take as No the group generated by Go, a reflection on a line parallel to b 
and reflections on the glide-reflection lines of Go and the midlines between any two 
consecutive glide reflection lines. Each fundamental domain of Go gets divided into 
8 fundamental domains of No, which are rectangles of sides j and k. See again 
Figure |l2](a) , where one of these fundamental subdomains has been shaded. 

Part (b) of Figure [l^ shows the extended Dirichlet region of a fundamental sub- 
domain, and part (c) the corresponding influence region. The shadowed rectangles 
in part (c) are the fundamental subdomains lying in the influence region but which 
cannot produce overlapping regions if the two orbits of Go are related by No, with 
the same argument as in part (a) of Figure |)| Hence, the interior of the thick poly- 
gon is the reduced influence region. We label the eight Go-cosets in No with the 
letters A, B, C, D, E, F, G and H where A is Go itself, and get for the orbits in each 
class the number of overlapping regions shown on the right of Figure ^(c). This 
number equals four in the cosets obtained from A by a translation in the direction 
of a (cosct C), reflection on a line parallel to b (cosets B and D), or glide-reflection 
with axis in the direction of a (coset G). It equals six in the coset E obtained by 
any of the other of reflections of Nq, and seven in the other two cosets F and G, 
obtained by order-two rotations or vertical translation. 

We now look at the two orthorhombic groups in Table ^ having a horizontal pg. 
They are depicted in Figure As usual, the picture shows the projection of a 
translational cell, which in this case consists of two fundamental domains of Go- In 
the group p2i2 2i we see that the horizontal planes contain orbits of types D, E 
and H, while in P^^IT we have of types D, F and G. 

Proposition 3.7. Dirichlet stereohedra for the groups P^f „ and P^^^r can- 
not have more than 38 facets, respectively. 

Proof: In -P^x^r, ^ ne band Zp contains two horizontal orbits of each of types D, 
F and G, which produce at most 2 x (4 + 7 + 4) = 30 neighbors. These, added to 
the six neighbors in the base horizontal orbit and the two vertical translates of P, 
provides the upper bound. 

For the group P^^^r we will use vertical planes, and Figure |TJ which shows 
the projection of a generic orbit of G to the reduced influence region computed 
above. 

Counting 2 points over each white rectangle and one over each grey rectangle 
gives 44 possible neighbors. 

Let us concentrate in the Go-orbit containing P, whose intersection with the 
reduced influence region consists of P and the eight points in grey rectangles, labeled 
X\, X 2 , Yi, Y 2 , Z\, Z 2 , Zj, and Z 4 in Figure O. X\, X 2 , Y\ and Y 2 are obtained 
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Figure 12. Extended Dirichlct region and influence region for 
the group pg 



• 






• 1/2 




• 






• 1/2 


o 






Ol/2 






01/2 


O 






Ol/2 


O 






O 






Ol/2 




• 1/2 


• 








• 1/2 


• 





n c a 



p2x 2lL 2lL 

a b c 



Figure 13. Orthorhombic groups with the horizontal planar 
group pg 
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FIGURE 14. Vertical planes in the group P^i^j- 



from P by translations of Go- The other four by glide reflections. In the planar 
Dirichlet tiling produced by this orbit we have: 

(i) P is a neighbor of the point X\ if and only if it is a neighbor of X%, because 
P being a neighbor of Q implies tP being a neighbor of tQ for any t £ Go. 
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(ii) P is a neighbor of the point Y\ if and only if it is a neighbor of Yi , for the 
same reason. 

(iii) The four points X\, X2, Y\ and Y2 cannot be all neighbors of P. This is a 
topological argument: If X\ and X2 are neighbors of P then Z\ is a neighbor 
of Z2 and Z 3 is a neighbor of Z4. There is no "room left" for Y\ or Y 2 being 
neighbors of P. 

Let us say that we are in the "X" case if Yi and Y2 are not neighbors of P and 
in the "Y" case if X\ and X2 are not neighbors of P. Statement (iii) above implies 
that we are always in one of the two cases (maybe in the two of them). 

The dashed lines in Figure [l4] represent four vertical planes R, S, X and Y. The 
number 44 implied counting 8 possible neighbors in each of them and 30 points in 
total in the four planes, since the two points counted over P are common to X and 
Y. We will see that the four planes can contain at most 26 neighbors of P in total, 
which finishes the proof. 



R and S can contain at most seven neighbors each by Theorem 1.4, since the 
planar subgroup of G with respect to the planes parallel to R, S and X is of type 
P9- 

X and Y contain P and hence each of them gives at most six neighbors. But we 
can be more precise: in the "X" case, X gives at most six neighbors and Y gives 
at most four neighbors not contained in X. In the "Y" case, Y gives at most six 
neighbors and X gives at most four neighbors not contained in Y. □ 



4. The group P6i22 

4.1. An upper bound using vertical planes. We now deal with a group G of 
type P6i22, whose horizontal subgroup Go is generated by two translations of equal 
length forming an angle of 60° . The whole group is generated by Go together with 

• a screw rotation of order 6 and of vertical axis and 

• any rotation of order 2 with axis parallel to one of the generating horizontal 
translations and intersecting the screw-rotation axis. 

See a graphical representation of the group in Figure |l5|(a). 




Figure 15. P6i22: (a) Neighbors in the interior of Z P (b) XZ view 



Proposition 4.1. Dirichlet stereohedra for the group P6±22 cannot have more than 
78 facets. 
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Proof: Remember that Corollary 1.6 gave an upper bound of 96 for the number of 
neighbors in this group. We hrst show a different way of deriving this same upper 
bound which will be more appropriate for our purposes here. 

Let I denote the minimal length of a horizontal translation in G, so that the 
horizontal group G generates a triangular lattice with equilateral triangles of side 
I. Each such triangle has height a = y/3l/2. Then, for any point P we have that 
VoiG a p(P) is an infinite prism over a regular hexagon of side 2a/3 = l/y/3. Since 
this happens for every point, the condition necessary and sufficient for Vor<3 p(P) D 
V° t G q(Q) ^ is than that Q lies in the prism over 2Vor g p(P), i.e., over the 
regular hexagon of side 4a/ 3 — 2Z/V3 centered at P and with sides orthogonal 
to the primitive translations of Go- This is the thick hexagon in Figure |l5|(b). 
Moreover, if Q is not in the same horizontal plane as P, a necessary condition for 
a Q £ GP being a neighbor of P in Vorcp(P) is that the projection of Q to the 
horizontal plane containing P lies strictly inside that hexagon. For the points in 
the same horizontal plane as P it is allowed, however, to lie in the boundary of the 
hexagon. Then, the possible neighbors of P in Vorcp(P) are: 

• The six points obtained from P by primitive horizontal translations, i.e., the 
mid-points of the edges of the thick hexagon in Figure [l5| 

• The two closer vertical translates of P. 

• The points which project to the interior of the thick hexagon and with vertical 
distance to P smaller than 1. There are 22 horizontal planes other than the 
one containing P and with distance to P smaller than 1. A simple density 
argument shows that at most four orbit points in each of the 22 planes lie 
inside the hexagon. More precisely, the density argument says that exactly 
four points at each horizontal plane project to the hexagon if one counts points 
projecting to the interior as 1, those projecting to facets of the hexagon as 
1/2 and those projecting to vertices as 1/3. 

Summing up, the above gives a bound of 22 x 4 + 6 + 2 = 96 neighbors. Our 
goal is to decrease the 22 x 4 part by considering vertical planes parallel to the 
XZ coordinate plane. The subgroup G\ consisting of elements of G which preserve 
those planes is of type p2, generated by a translation of length I in direction X, 
a translation of length 1 in direction Z and any order 2 rotation contained in Gi 
(whose axis will be parallel to Y and intersect some screw rotation axes of G). 

The orbit GP decomposes into an infinite family of orbits of the group Gi lying 
in different X Z-planes. Ten of them are marked in Figure |l^(b) with the letters A, 
B, G, D, E, A' , B', C', D' and E 1 . A is the one containing the base point P. In each 
pair A- A', B-B' , C-C', D-D' and E-E' one of the G x orbits is obtained from the 
other by a screw rotation p of order 2 in the vertical axis {X = —1/4, Y = —a/2}. 
More importantly, when projected to any of the XZ-planes, p becomes an element 
of the normalizer of G\ which exchanges the coloring of the rectangular tiling used 



in Lemma 2.3. Hence, by that Lemma, 



Claim 1: In A' we can have at most 4 neighbors of P and in any of the other 
pairs B-B', C-C', D-D' and E-E' one of the two planes can provide at most 7 
neighbors and the other one at most 4 neighbors (although a priori we do not know 
which plane is which). 

Let us now show that the bound of 96 obtained above over-counted the points in 
these ten planes by at least 18. For this we use Figure |l^, which essentially coincides 
with Figure |l5|(b) except that we have removed the numbers showing the height of 
different points and included instead the subdivision of the plane into fundamental 
subdomains of the horizontal group Go- The darker fundamental subdomain, which 
we call Do, is the one containing the base point P. 
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The lighter shaded polygonal region of Figure |l6| consists of the fundamental 
subdomains D with the property that for any choice of P in Dq the points of GP 
lying over D are in the interior of the hexagonal influence region centered at P. 
In other words, the fundamental subdomains D for which the upper bound of 96 
includes two points above D for any choice of P. 




Figure 16. Vertical planes in the group P6i22 

This implies that the 96 points of the previous bound include at least 6 points in 
the planes D' , C and E, at least 8 points in the planes B' , A' , A, B, C and D and 
at least 5 points in the plane E' . The number for A include the special counting of 
points obtained from P by translations of G. The number for E 1 takes into account 
a fundamental subdomain whose orbit points lie in a facet of the hexagonal prism 
and which have been counted as 1/2 in the "96 points". According to Claim 1, 
we have over-counted at least 5, 3 and 3, respectively, in the pairs of planes B-B', 
C-C and D-D'. In the pair E-E' we have over-counted at least 1, since one of 
them produces at most 4 neighbors and we have counted more than 4 in each of 
them. Finally, we have over-counted 4 points in the plane A' and 2 in the plane 
A. □ 

4.2. Dirichlet stereohedra with many facets. Our constructions of Dirichlet 
stereohedra with many facets are based on the following result. The case k = 8 
appeared in . The general case has been independently proved by J. Erickson as 
Lemma 1 in [p|: 

Lemma 4.2. Let g be the screw rotation of order k around the Z coordinate axis, 
with translation of length I, i.e. g{x,y, z) = (x cos(2n / k) ~ y sin(2ir / k) , x sm(2ir / k) + 
y cos(27r/fc), z + I). Let (g) be the (infinite cyclic) group generated by g. Then, any 
point P not on the axis of g has as neighbors in Vor^p(P) the 2k points g l P for 
i€{-k,...,k}\{0}. 

Proof: For any P £ R 3 (but not on the Z-axis) the convex hull of (g)P is an 
infinite prism over a regular A-gon, one of the edges of the prism containing the 
points g nk P, for n £ Z. In particular, g k P and g~ k (P) must be neighbors of P. We 
will prove that the points g z P for i £ {— k + 1, . . . , k — 1} \ {0} are also neighbors 
of P. 

For any P £ R 3 the orbit (g)P is contained in an helicoidal curve H r a i '. — 
{(rcos(t — a),rsin(< — a), j^t) £ R 3 : t £ R}. There is no loss of generality in 
fixing I — 2tt and a — 0, hence using the notation H r for H r 2 -n- We denote 
Pt = (r cos(i), r sin(i), t) the point of H r for a certain parameter t, so that g l Pt = 
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Pt+2-Ki/k- We will prove that for any two points Pt and Pp whose angular distance 
\t — 1'\ is less than 2tt, there is a sphere tangent to H r at these two points and not 
containing any other point of H r . This implies the lemma. 

To prove our claim, by symmetry considerations we can further assume that t' = 
—t. Then, the two points in the claim are of the form P to = (r cos(to), r sin (to), to) 
and P-t — ( r cos(to), — r sin (to), — to) for a certain to € (0, tt). The sphere tangent 
to H r at P to and P_ to will have center at a point O — (a, 0, 0) which satisfies that 
the vectors 

OP to = (rcos(to) - a,rsin(t ),to) 

and 

(dP t /dt) t= t = (r sin(t ), -rcos(to), 1) 
are orthogonal. The equation 

(r cos(t ) — a)r sin (to) — r 2 cos(to) sin(t ) + to = 
gives the solution 

-to 



rsin(to) 

We now prove that for any point Pt other than P to or P~t a the distance from 
O to P t is strictly bigger than the distance from O to Pt . It suffices to consider 
t € [ — 7r, 7r] because d{0,P t ±2Tr) > d(0,P t ) for any t G [— 7r, 7r]. Moreover, since 
d(0, Pt) — d(0, P-t) we restrict our attention to t e [0, tt]. We define the function 

/(t) := d(0, P t f - d(0, P to f = ^(cos(t)-cos(to)) + f _ 

sm(to) 

whose first and second derivatives are 

i . ° . sin(t), , 
sm(t ) sm(to) 



/'(*) = 2* - sin(t), /"(t) = 2 - 2^- cos(t). 



We have 

(i) f(t Q ) = /'(to) = (as expected). 

(ii) /"(to) > 0. For this, observe that /"(to) > is equivalent to sin(t ) — 
to cos(to) > 0, which holds because the function g(t) := sin(t) — t cos(t) is zero 
at the origin and its derivative g'(t) = tsin(t) is strictly positive in (0,7r). 

(iii) /' is injective in (0,7r). Indeed, if /'(ti) = f'fa) then sin(ti)/ti = sin(t 2 )/t2 
and, in particular, the derivative of the function h(t) := sin(t)/t must have 
a zero between ti and t 2 . But h'(t) = (tcos(t) — sin(t))/t 2 = would imply 
t = tan(t), which does not happen in the interval (0, tt). 

With this we prove that / achieves its unique minimum in the interval [0, tt] 
at f = to, as follows: Claim (iii) implies that /" is either always non- negative or 
always non-positive. By (ii) it is always non-negative. Claims (i) and (ii) imply 
that / has a local minimum at to which, by the previous observation, is the unique 
global minimum. □ 

Let G be a crystallographic group of type P6i22. Its two metric parameters 
(which define G modulo conjugation by an isometry) are the lengths I and a of its 
minimal translations in directions parallel and perpendicular, respectively, to the 
order-6 screw-rotation axes. As usual, we assume that the screw rotation axes are 
vertical, that is to say, parallel to the third coordinate axis. Also, assume that one 
of the minimal horizontal translations is parallel to the X-axis and that the X-axis 
itself is one of the horizontal order-2 rotation axes in G. This completely specifies 
the group G (and agrees with Figures [If] and |l6| ). 

Let g be the screw rotation of order twelve obtained substituting k = 12 in 
Lemma |4.2|. Observe that g 2 € G. Let p be the order-2 rotation (x, y, z) ^ 
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(x, —y, —z) which is in G by hypothesis. Suppose now that our base point P has 
coordinates (rcos(7r/12), r sin(7r/12), Z/24), so that p(P) = g^ 1 (P). Under these 
assumptions we have that the orbits of P under the groups (g) and (p, g 2 ) C G 
coincide. By Lemma [h^, the Voronoi diagram of this part of GP alone produces 
24 neighbors of P. Now, if we fix the parameters I and r in the above description 
and make a tend to infinity, the neighbors of P in the Voronoi diagram of (p, g 2 )P 
will keep being neighbors in GP. (Observe that GP is obtained as the Minkowski 
sum of {p,g 2 }P and a triangular grid of side a in a horizontal plane). Moreover, 
since the regions in the Voronoi diagram of (p, g 2 )P are unbounded, new neighbors 
are guaranteed to appear. Hence, the above construction is guaranteed to produce 
at least 25 neighbors. 

We could try to argue geometrically how many "new" neighbors have to appear 
in the limit of a going to infinity. For example, it is relatively easy to show that 
this number does not depend on the choice of / and r. However, it seems easier to 
compute that number experimentally: 

Example 4.3. Taking I = 12, r = 1 and a = 100 in the above setting (which gives 
P = (cos(7r/12), sin(7r/12), 1/2)) one gets 31 neighbors. The number and identity 
of the neighbors is (experimentally) stable under increasing the value of a. 

Even more, playing with different possibilities for the parameters we have found 
that: 

Example 4.4. The metric parameters I — 100, a — 950 and the base point P = 
(1, tan(7r/12), 4) produce a Dirichlet stereohedron with 32 neighbors. This stere- 
ohedron is very unstable. For example, changing the last coordinate of P to be 
100/24 = 4.166, which would match exactly the above setting with r — 1/ cos(7r/12), 
only 30 neighbors are obtained (and their identity changes drastically). Also, the 
parameter a is not big enough for the number of neighbors to be the same for bigger 
values of a. 

Figure [l?] describes the points of the orbit which produce neighbors in each of 
the two examples. It shows the same decomposition of the plane into fundamental 
subdomains which appears in Figure |l^. The grey subdomain is the one containing 
P. The points of the helicoidal orbit (p,g 2 )P are those projecting to the regular 
hexagon. A plus (resp. a minus) in a subdomain means that the point of GP 
projecting to that subdomain and lying in the upper (resp. lower) half of the band 
Zp is a neighbor of P. In the first example all the 24 possible neighbors within the 
hexagon appear, as predicted. In the second example only 16 of them appear, but 
this is compensated by more neighbors out of the hexagon. 




(a) 



(b) 



Figure 17. Neighbors of P obtained in Examples [13 
(parts (a) and (b) respectively, for a group P6i22. 



and 



4.4 
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The computations were made using a Maple program which first generates a 
sufficiently large number of points S of the orbit GP to guarantee that all the 
neighbors of P lie in S and then checks for each of them whether it is actually a 
neighbor or not. For the first step, it would be enough for example to use the two 
points in Zp in each of the subdomains which appear in Figure [l(| For the second 
step, we express being a neighbor of P as feasibility of a certain linear program. 

The same ideas can be applied to groups having an order-4 screw rotation around 
a vertical axis and an order-2 rotation in a horizontal axis, with the two axes in- 
tersecting one another. One has to use Lemma 4.2 with k = 8 and, hence, can 



only guarantee to obtain more than 16 neighbors. This is the way we constructed 
stereohedra with 18 facets for the group Ig^^i hence showing that this is ex- 
actly the highest possible number of facets of Dirichlet stereohedra for groups with 
reflections (see ||, Example 2.9]). 

Specially good is the group 74i22, whose graphical representation appears in 
part (a) of Figure [l8|. It has screw rotations or order-4 in the two versions 'dextro' 
and 'levo', and allows a point to be considered as lying in two different 'helices'. 

Example 4.5. The following base point and metric parameters for a tetragonal 
group 74i22 produce Dirichlet stereohedra with 29 facets: 

Minimal length of horizontal translation = 4 
Minimal length of vertical translation = 1 
Base point P = (1, -|, j^), in the coordinate system of Figure [7^ (a) 

Part (b) of Figure [l8| shows the 29 orbit points producing facets of the stereohe- 
dron, with the same conventions of the previous examples. The zero in one of the 
fundamental subdomains indicates that the corresponding neighbor is at the same 
height as the base point P. 
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Figure 18. Neighbors of P obtained in Example 4.5 for a group 74i22 
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